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THEORETICAL AND EXPERIMENTAL POWER FROM LARfiE 

HORIZONTAL-AXIS WIND TURBINES* 

Larry A. Vlterna and David C. Janetzke 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 


SUMMARY 

A method for calculating the output power from large horizontal-axis wind 
turbines is presented. Modifications to the airfoil characteristics and the 
momentum portion of classical blade element-momentum theory are given that 
Improve correlation with measured data. Improvement is particularly evident 
at low tip speed ratios where aerodynamic stall can occur as the blade experi- 
ences high angles of attack. 

Output power calculated using the modified theory is compared with measured 
data for several large wind turbines. These wind turbines range In size from 
the DOE /NASA 100 kW Mod-0 (38 m rotor diameter) to the 2000 kW Mod-1 (61 m 
rotor diameter). The calculated results are In good agreement with measured 
data from these machines. 


INTRODUCTION 

To correctly assess the cost of energy of the next generation of wind tur- 
bines, a validated method of predicting aerodynamic performance Is required. 
Several methods have been available (refs. 1 and 2), however, the results from 
which have not been sufficiently verified with experimental data from a variety 
of operating wind turbines. 

Classical blade element-momentum theory has been found to give adequate 
results under normal operating conditions near the design tip-speed ratio 
(ratio of blade tip speed to wind speed). At both high and low tip-speed 
ratios, however, the wind turbine may operate in flow conditions which are not 
easily analyzed by classical theory. Recent tests on the NASA Mod-0 100 kW 
wind turbine, for example, indicate that classical theory is inadequate at low 
tip-speed ratios when the airfoils are at high angles of attack. This problem 
can be particularly important in the calculation of fixed-pitch and tip-control 
rotor performance. Since the maximum power produced by a fixed pitch rotor is 
a critical design parameter, efforts are being made to improve theoretical 
analysis of this operating condition. At high tip-speed ratios, the rotor may 
operate in the vortex ring state in which a vortex ring is formed in the rotor 
plane. Under this condition an eddy motion occurs, such as behind a bluff 
body, and the overall thrust force On the rotor can deviate significantly 


*8ased on a presentation given at the Fifth Biennial Conference and Workshop 
on Wind Energy, October 5-7, 1981. SERI/CP-635-1340, DOE Publication 
CONF-811043, vol. 2. 


from valuos given by momentum theory (ref. 3). Furthermore* from 
^ viewpoint, analysis of this state using momentum theory can cause 
ilUticult convergence problems and result in unrealistic values in the calcula- 
tion of the iruluceil flow velocities, 

Experimental data from two Mod-0 rotor configurations as well as from the 
Mod-OA and Mod-1 wind turbines were analyzed. Empirical corrections to both 
the airfoil characteristics and momentum theory are presented which improve the 
correlation with experimental results. Improvement is most noticeable in the 
power output versus wind speed at low tip speed ratios. 


CLASSICAL THEORY 

Blade element-momentum theory used in the PROP Code (ref. 1), as well as 
others, divides the blade Into small spanwise elements. These elements are 
considered to act as airfoil segments in two-dimensional flow fields, each at 
a particular angle of attack. The effects of tower shadow and nonuniform in- 
flow have been found to be small and are generally neglected. From the geome- 
try, the rotational velocity component, the wind component, and the ''induced" 
axial and rotational components, the local angle of attack is calculated. The 
lift and drag forces on the elements are then determined from two-dimensional 
(infinite aspect ratio) wind tunnel data at the local angle of attack (ref. 4). 
Comparison of cambered and uncambered data indicates little difference in air- 
foil performance beyond stall; therefore, data for synynetric sections may be 
used for this region (ref, 5). The "induced effects" of the wind turbine on 
the flow are determined in an iterative procedure until momentum considerations 
are satisfied. The end effects of the finite-length wind turbine blade are in- 
cluded by using a tip loss model. There are a variety of these tip-loss models, 
the simplest being to reduce the lift coefficient to zero tor approximately 3 
percent of the radius near the tips. A more complete description of blade 
element-momentum theory is contained in reference 1. 


LOW TIP SPEED RATIO OPERATION 

Tests were run on the NASA Mod-0 100 kW wind turbine to investigate the 
performance characteristics of fixed-pitch rotors at low tip-speed ratios (high 
winds). Under this Condition, the airfoils operate at high angles of attack 
and aerodynamic stall can occur. 

Test Configurations and Operating Conditions 

Two rotors were used in the Mod-0 tests, both of which were 38 m in diam- 
eter. Details of the blade plz.tform, twist, etc. are given in tables 1(a) and 
(b) and figures 1(a) and (b)i The significant differences between the two 
blades are that the aluminum blades have a variable thickness-to-ciiord ratio 
(Wc) and 34 of nonlinear twist, while the steel spar blades have a constant 
t/c and no twist. The thickness and twist distributions for the aluminum 
blades are given in figures 2(a) and (b). 


The Mod-0 wind turbine was operated with the highly twisted aluminum blades 
at nominal rotor speeds of 20 and 2b rpin. Due to slip of the fluid coupling in 
the drive tram, actual rotor speeds were 21.0 and 27.4 rpm at the maximum powers 
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of 50 kW and 105 kW respectively. The blades remained fixed in pitch throughout 
these tests. Alternator output power and reference windspeed were recorded on 
magnetic tape and analyzed using the method of bins (ref. 6), A reference wind- 
speed is measured by one of five hub height anemometers located in a semicir- 
cular array of 60 m radius around the wind turbine. The anemometer, which is 
most directly upwind, is chosen as the reference since it is believed to give 
windspeeds which are representative of the average freestream windspeed near the 
rotor. 

Further tests were conducted using the steel spar blades with no twist. The 
rotor speed was 32 rpm. The inboard 70 percent of these blades remain fixed in 
pitch and thus experience high angles of attack at high windspeeds. As the wind 
speed increases, the wind turbine power output increases until the generator 
rating of 100 kW is reached. When wind speed increases further, it is necessary 
to pitch the outboard blade section so that the generator rating will not be ex- 
ceeded. Since data from the anemometer array was not available at this time, a 
nacelle wind speed was used. The nacelle wind speed is measured by an anemom- 
eter located on the wind turbine nacelle. 

Baseline Performance Analysis Model 

The Mod-0 rotors were modelled with the PROP Code. The aerodynamic char- 
acteristics used for the aluminum blade are given in figure 3 for a NACA 23018 
''half-rough'' airfoil. The characteristics shown in figure 3 are not measured 
data but instead were derived from a number of sources using some engineering 
judgment. The designation "half-rough" denotes aerodynamic data which are an 
average of NACA smooth and rough data given in reference 4, This roughness ef- 
fect was included to account for manufacturing imperfections and for the accumu- 
lation of dirt as the airfoil is exposed to the environment. 

Since the output of the PROP Code is rotor power with no drive-train los- 
ses, the following model based on experimental data was used to calculate the 
electrical power: 

P3 « -1.9 + 0.82 P2 (1) 

in which P3 is the electrical power, (kW) and Pg is the rotor power, (kW). 

Comparison of Measured and Predicted Performance (Baseline Model) 

Figure 4 shows the predicted and measured power versus windspeed for the 
Mod-0 aluminum blades. The correlation between measured data and analysis 
using standard two-dimensional airfoil data is not very good. The obvious 
deficiencies of the predicted results are the less-than-measured maximum power 
and the rapid decrease in power after stall. The most interesting character- 
istic of the data is the relatively constant output power at high windspeeds 
for both rotor speeds. The leveling off of the power output occurs for wind 
speeds at which stall occurs over most of the blade. Maximum power is greater 
for the higher rotor speed since higher wind speed is required to achieve the 
same stall angle of attack. Results obtained for the Mod-0 steel spar blade 
are shown in figure 5. As shown in the data, the outboard section pitching 
continues to increase as wind speed increases. This data indicates a contin- 
uing production of torque by the inboard fixed-pitch portion of the blade, 
even though that section is stalled and producing negative power according to 
the model. 
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Improved Performance Analysis Model 

In an attempt to improve the correlation with experimental data, smooth 
airfoil characteristics were used. The fact that the baseline analysis under- 
predicted the actual performance suggested that the use of smooth airfoil data 
might improve correlation. By incorporating smooth airfoil characteristics, an 
increase in the maximum power resulted from the increased maximum lift and low- 
er drag. Thus, it is recommended that smooth airfoil data be used to analyze 
wind turbine rotor performance near stall. 

The use of smooth airfoil data, however, did not improve the rapid decrease in 
power after stall. Investigation of the airfoil characteristics reveals the 
reason for the predicted negative power at very high winds. Figure 6 shows an 
airfoil element operating with its chord line parallel to the plane of >*ota- 
tion, such as on the inboard portion of the steel spar rotor. The resultant 
wind velocity acts at an angle of attack, o, with respect to the plane of ro- 


tation. The coefficients of lift. Cl, and drag, Cn, forces operating on 
the element can be resolved Into a coefficient of torgue force, Cq , which 
acts in the plane of rotation. This coefficient is given by 

Cq « sin o - Cp cos o (2) 

Using the baseline two dimensional airfoil data from figure 3 in equation 
(2) produces the torque coefficient characteristic shown in figure 7. This 
figurG shows that negdtivG torguG can be expected for angles of attack between 
stall and 45*. 

To improve the correlation between the measured and calculated performance 
after stall, an idealized stall model was used. This model was developed by 
deriving the airfoil characteristics in the stalled region that would produce 
performance similar to the measured data. The conditions placed on the model 
were threefold: (1) continuity of the stalled and unstalled characteristics 

at the stall angle, (2) constant power (or torque coefficient) after stall, 
and (3) the correct drag coefficient at 90 angle of attack for the given 
blade aspect ratio. The derivation of this model follows. 

The torque force on a blade element is proportional to the coefficient of 
torque force and the square of the resultant velocity, or mathematically, 

Q ' 

For a constant rotor speed we can divide by the constant V^, which yields 

g -V Cp ^ (4 

0 

but from figure 6, cos a = 


but from figure 6, 


and thus 


cos o 
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Substitutinu Equation, (2) gives 

sin a Cp 

g n, — j- — - — -~- 

cos'^o COS a 


(6) 


Two particularly convenient forms for the lift and drag characteristics 
after stall that satisfy the conditions stated above are given by 


AjSin 2 a + Ag 



(7) 


and 

Cp » B^Sin^a + BgCOS a ^8) 

Substituting these into Equation (6) yields 

Q <v (2A^ - B^)sin o tan o + (Ag - Bg) (9) 

Since the measured torque after stall is independent of wind speed it must 
also be independent of angle of attack. Thus, knowing that the derivative with 
respect to a must be zero, by inspection 



At an angle of attack of 90*, equation (8) lives 

h “ ^DMAX 


where Cp|v|/\x Is the drag coefficient at a = 90*. 
For a finite aspect ratio blade 


'"DMAX “ * 0.18« 


( 10 ) 

( 11 ) 


( 12 ) 


for 50 based on experimental data from reference 7. 
Thus 


ft ^DMAX 

A « — 

Rearranging equation (7) and substituting equation (13) yields 

sin a 


A., 


(C^ — ^pjvy^X ® cos a) 


T 

cos a 


(13) 


(14) 
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Slmlldfly equation (8) and (ll) 9 ^ve 

Cn - C 


8, 


0 ’ ‘'DMAX 


sin^a 


cos a 


(15) 


For continuity with the below stall airfoil data, Ag and B? 
solved for the stall angle conditions, and thus equations (14) and (15) become 


^2 “ ■ ^OMAX ®s 


sin 0 . 

cos a-) 


"s' “T 
^ cos a, 


(16) 


and 




^DMAX 
cos Be 


(17) 


In which the subscript s denotes the value of the constant at stall. 

The resulting airfoil characteristics for a NACA 23018 
in figures 8 and 9 along with the standard two-dimensional airfoil data from 

figures 3 and 7. 

Comparison of Measured and Predicted Performance (Improved Model) 

With these modifications to the airfoil characteristics, 

Mod-0 Steel spar blade as shown in figure U. Since t"®sej:wo oiaoes are 
sionif icantly different in planform and twist, it is expected that this 
approach win produce satisfactory results for other configurations. 

HIGH TIP SPEED RATIO OPERATION 

At hiah tip-speed ratios (low winds), a wind turbine may operate 
vortex ring sUter Under this condition, the flow "»l6l developed by blade 
element-momentum theory breaks down. The induced axial interference fkitor, 
a ^SveraertrvalUM greater than 0.6. Since the flow far downstream would 
have twicl the induced velocity, this suggests J ”1; "fJ'J Jhe «er!' 

It ';of?hrls;^rfliJieit"Si)fe?s®g^^^^^^^ ?L'';ii;erJal'c;uud using 

momentum theory as shown in figure 12. The 8) 

calculation of the local thrust coefficients. 

that the following linear algebraic expression for the local thrust coeffl 
cient may be used: 

C| = 0.5775 + 0.96a a < 0.38 (15) 

The use of this empirical equation can eliminate convergence P^^oblems of 
the induced axial interference factor at high tip-speed ratios. Figure 13 
shows a graphical representation of how the value of the interference factor 
is determined such that both momentum and blade element theories are satis- 
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fled. For a tip-speed ratio of 8» the two theories converge for a value of 
"a" of 0.37. At a tip-speed ratio of 12, however, no Intersection occurs be- 
tween the curves representing momentum and blade element theories. The use of 
the empirical curve, though, gives a converged solution at an "a" factor of 
0.58. 

Figure 14 shows measured power In low winds (high tip speed ratios) for 
the Mod-1 wind turbine described In table II. The Mod-1 Is presented In this 
example since It Is more likely to operate In the vortex ring state due to Its 
high rotor solidity and high tip speed ratio. Theoretical power output using 
four methods Is shown for comparison. The solid curve shows the results using 
the standard PROP code where the Interference factor Is limited to 1. This 
model tends to fall near the lower bound of the measured data. Limiting the 
Interference factor to 0.5 has also been used since It prevents theoretical 
flow reversal In the far wake. However, as shown in figure 14, that model 
predicts power levels near the upper bound of the data. Limiting the Inter- 
ference factor to 0.7 or using the empirical curve gives results which more 
nearly represent the average of the data. For simplicity the Interference 
factor was limited to 0.7 In the calculations of the next section. 


FURTHER COMPARISON OF THEORY WITH FIELD MACHINE DATA 

For verification of the PROP code, comparisons were made with measured 
data from the Mod-OA and Mod-1 wind turbines. The Mod-OA and Mod-1 wind tur- 
bines are field machines which were designed to operate In an unattended mode. 
Engineering data Is automatically recorded on magnetic tape. Ten-mln averages 
of the power and wind speed are calculated using an energy method described In 
reference 9. 

Figure 15 shows the measured and theoretical power output versus wind 
speed for the Mod-OA blades. The blade planform characteristics are given in 
table 1(a) and figure 2. With the exception of at the rated wind speed (9.5 
m/s), the measured data fall both above and below the theoretical curve. Near 
rated wind speed the data fall below the theoretical curve since the power Is 
an average of times both below and at rated power. The mean power output den- 
sity (power/rotor swept area) for the entire data set as 90 ^/wr. The theo- 
retical mean power output density, however, was 89 W/m*- for a mean deviation 
of +1 percent. Figure 16 shows the measured and theoretical power output for 
the Mod-1 60 m diameter steel blades described In table II. The measured data 
from this machine generally fall above the theoretical power output curve. The 
mean power output density for this data set was 250 W/nr. The theoretical 
mean power output density, however, was 238 W/m^ for a deviation of +5 percent. 

CONCLUSIONS 

A method for calculating the output power from large horizontal axis wind 
turbines has been developed. This method Is based on blade element-momentum 
theory with empirical modifications to the airfoil data and momentum theory. 
Smooth airfoil data Is used for the unstalled operating region. Post-stall 
airfoil characteristics were derived to improve correlation with measured data, 
particularly at low tip-speed ratios. At high tip-speed ratios, empirical 
equations are used Instead of momentum theory to calculate the values of the 
local thrust coefficients. 
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Output power calculated using the Improved method were compared with data 
from several large wind turbines ranging In size from the 100 kW OOt/NASA 
Mod-0 (3B-III rotor diam) to the 2000 kW Mod-l (61-m rotor diaiii). The theoret- 
ical mean power densities are within 5 percent of the measured values for the 
cases exaiii1rie<l. 


a 

A 

Ai. A2 
Bi. B2 
B 
c 

Cq 

^Dmax 

Cl 

Cq 

Ct 

P2 

P3 

Q 

r 

R 

Vr 

Vo 

Vl 

V 

a 

a 

s 


NOMENCLATURE 

axial Interference factor, a o (Vo-Vi)/Vo 

aspect ratio of blade, R/c, based chord length at 75 percent 
radius 

constants In lift coefficient equation after stall 

constants In drag coefficient equation after stall 

number of blades 

chord length, m 

drag force coefficient 

drag force coefficient at 90* angle of attack 

lift force coefficient 

torque force coefficient 

thrust force coefficient 

rotor output power, kW 

electrical output power, kW 

torque force, N 

local radius, m 

radius of rotor blade, m 

resultant velocity vector, m/s 

free-stream wind speed, m/s 

wind velocity at rotor plane, m/s 

velocity due to rotation, m/s 

angle of attack, deg 

blade twist angle, deg 

value of airfoil characteristic at stall 
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TAliUE 1(a), " MUP"0/-0A AUUMINUM B1.API; CHARACTERISTICS 


Rotor clla., in(ft) 

Number of blades 

Root cutout, percent span 

Pitch control (U^A) . . , 

Ulade^^itch 7b% span, deg (0/UA +90 ■ Feather) 

Taper ’ ,* * 

Twist, (leg 

Solidity , , 

Precone, deg (O/OA) 

Rotor shaft tilt, deg (O/OA) 


. 30,5 (li^5.4) 

...... i 

f Ixed/varlable 
. . . - 2 . 8 / 0.0 
NACA 230 series 
.... linear 
34 (nonlinear) 
.... 0.031 
. . . 0.0/7.0 
. . . 8. 5/0.0 


TABLE 1(b). - MOO-0 STEEL SPAR BLADE CHARACTERISTICS 


Rotor dia., m(ft) 

Number of blades ....... 

Root cutout, percent span . . . 

Pitch control 

Blade pitch Inboard section, deg 

Airfoil . 

Taper , 

Solidity 

Twist, deg 

Precone, deg 

Rotor shaft tilt, deg 


38.4 (126.0) 

2 



30 percent span tip control 

0 

NACA 23024 

linear 

0.033 

0 

0 

8.5 


TABLE II. - MOO-1 BLADE CHARACTERISTICS 


Rotor dia., m(ft) 

Number of Blades 

Root cutout, percent span 

Pitch control 

Blade pitch 75% span, deg (+90 = Feather) 

Airfoil 

Taper 

Twist, deg 

Solidity 

Precone, deg 

Rotor shaft tilt, deg 


. 61 (201.6) 

2 

10 

. . variable 

NACA 44 series 
. . . linear 
. 11 (linear) 
.... 0.043 

• . • • . 9 

• . . • . 0 
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Figure 10* - Com]>erlson K (t^tasuredand 
calculated performance for the Mod- 0 
aluminum blades* 



Figure 12. * Comparison of theoretical and measured 
rotor thrust force coefficients. 
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Figure 13. * Theoretical and empirical rotor thrust force 
coefficients. 
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Figure R - bperlmenkl end theoretical output power from the Wod* 1 
wind energy system. 




